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Bulk-mediated surface diffusion on a cylinder: propagators and crossovers 
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We consider the effective surface motion of a particle that freely diffuses in the bulk and inter- 
mittently binds to that surface. From an exact approach we derive various regimes of the effective 
surface motion characterized by physical rates for binding/unbinding and the bulk diffusivity. We 
obtain a transient regime of superdiffusion and, in particular, a saturation regime characteristic for 
the cylindrical geometry. This saturation, however, in a finite system is not terminal but eventually 
turns over to normal surface diffusion. The first passage behavior of particles to the cylinder surface 
is derived. Consequences for actual systems are discussed. 
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Bulk mediated surface diffusion (BMSD) defines the ef- 
fective surface motion of a particle on a reactive surface 
that intermittently unbinds and diffuses in the adjacent 
bulk before rebinding (Fig. [T|). BMSD was revealed by 
NMR in porous glasses [ij and has relevance to numer- 
ous technological applications Q. The particular case 
of BMSD on a cylindrical surface is of importance for 
facilitated diffusion in gene regulation 0, H[ , the net mo- 
tion of motor proteins along cytoskeletal filaments 
the transient binding of chemicals to nanotubes 0, or 
the exchange behavior between cell surface and surround- 
ing bulk of rod-shaped bacteria (bacilli) and their linear 
arrangements [3] to name but a few examples. 

BMSD was previously investigated for a planar sur- 
face in terms of scaling argumen ts B , §], master equa- 
tion schemes and simulations [10( . More recently the 
first passage problem between particle unbinding and re- 
binding for a free cylindrical surface was considered [Til ] . 
Here we establish an exact treatment of BMSD for a re- 
active cylindrical surface deriving explicit expressions for 
the surface occupation, the effective mean squared dis- 
placement (MSD) on the surface, and the returning time 
distribution from the bulk. Different regimes ranging 
from transient superdiffusion to terminal normal diffu- 
sion emerge naturally from the physical timescales en- 
tering our description. We derive previously unknown 
regimes characteristic of the cylindrical geometry, most 
remarkably the existence of extremely long jumps as well 
as a saturation regime of the surface MSD that will be 
crucial to fully appreciate effective surface motion medi- 
ated by bulk diffusion and its experimental investigation. 

We consider a cylinder with radius a centered along 
the symmetry axis of an outer cylinder of radius b. As 
we are interested in the effective displacement along the 
cylinder in z direction (Fig. [1]) we assume rotational sym- 
metry around the z axis. The concentration of particles 
(density function, DF) in the bulk is denoted by C(r, z, t) 




Figure 1: A particle diffuses in the bulk (full lines) and in- 
termittently binds to a surface on which it can also diffuse 
(broken lines). This produces an effective surface motion. 



([C] = cm 3 ) while n(z, t) ([n] = cm *) refers to the par- 
ticle DF per length on the inner cylinder surface. The 
particle numbers are Nb(t) — 2tt f^rdr C(r, z,t)dz 
in the bulk and N s (t) — J x oo n(z,t)dz on the cylinder 
surface. With surface and bulk diffusivities D s and D b 
the two densities evolve according to the diffusion equa- 
tions (a derivation from a discrete model will be pre- 
sented elsewhere LIS 
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in axial direction on the cylinder and 
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C(r,z,t) (2) 



in the bulk using cylindrical coordinates. Note the cou- 
pling term between surface and bulk DFs through the 
flux term in Eq. ([T]) . This term is positive when the par- 
ticle concentration is higher above the cylinder surface 



2 



and negative otherwise. We equip the diffusion equations 
with the following boundary conditions: 



lim C(r, z, t) = fin(z, t), 



, dC 
and — — 

or 



= 0. 



(3) 



Thus right above the cylinder surface the bulk concentra- 
tion is defined by the surface density where the coupling 
constant [i = l/(2nakbT g) involves the binding rate kb 
and the mean unbinding time t q. The second relation 
defines a reflecting condition at r — b. The initial condi- 
tion corresponds to a sharp concentration at the surface: 



Umnfz, t) = NqS(z) and lim C(r, z, t) = 0. 
t->o t->o 



(4) 



Averaging over the diffusion equations we readily see that 
N s (t)+N b (t) = N . 

Fourier-Laplace transforming the DFs according to 

/oo poo 
dze lkz / dte- ut C(r,z,t) (5) 
-oo JO 

and analogously for n(k,u) (we use the explicit depen- 
dence on the Fourier and Laplace variables to denote the 
transforms), the surface DF yields in the form 



n(k, u) = 



Nn 



k 2 D s + ftgAi/A 



(6) 



where q = y/k 2 + u/D b and k = 2nafiD b — D b /{k b T Q s)- 
We thus find in the denominator of the usual diffusion 
propagator a correction proportional to k stemming from 
the bulk exchange. For the bulk DF we obtain 



C(r, k, u) = [iNq 



Kx(qb)I (gr) + h{qb)K {qr) 



A[u + D s k 2 + KqAx/A] 
where the A-factors are defined as follows: 

Ai = KifaaMqb) - hiq^K^qb), 
A = Ioiq^Ktiqb) + h(qb)K (qa) 



(7) 



(8) 



in terms of the modified Bessel functions. 

We note that the coupling strength k is a measure for 
the efficiency of the bulk-surface exchange. At k = the 
number of surface particles N s (t) remains constant. We 
define the coupling time scale 



*„ = D b /n 2 = [k b T oS } 2 /D b 



(9) 



that tends to infinity for vanishing coupling and to zero 
for strong coupling (k — > oo). From the diffusion behav- 
ior we extract two additional timescales, namely 



t a = a 2 /D b and t b = b 2 /D b 



(10) 



For times shorter than t a the particle does not yet sense 
the curvature of the cylinder surface while for times 



longer than t b it feels the confinement by the outer cylin- 
der. In what follows we are mainly interested in the 
regime of strong coupling, t K -C t a -C t b but will also re- 
port an almost ballistic behavior under weaker coupling. 
The remaining cases will be discussed elsewhere |12j . 

The number of adsorbed particles follows from n(k, u) 
by taking k — > 0. For short times t <C t K we find that 
N s (t) w N remains constant as it should by definition 
of the coupling time. At longer times t K -C t -C t a the 
behavior changes to N s (t) sa Noy/t K /[Trt] ~ t~ x / 2 . For 
even longer times t a <C t -C h we have a faster decay 
N s (t) » \Noy/t K t a /t inversely proportional to t. Finally 
at very long times t > tb the dynamics equilibrates with 
respect to the radial diffusion and N s (t) ~ 2NQ^Jt K t a /tb. 

Surface diffusion. The MSD on the cylinder surface 
can be obtained from the characteristic function through 



{z\u)) 



-No 1 



d 2 n(k, u) 



dk 2 



(11) 



k=0 



where we divide by Nq to obtain an effective one particle 
displacement. This quantity is not corrected for particles 
leaving the surface in contrast to the normalized displace- 
ment {z 2 (u)) n = N (z 2 (u))/N s (t). For the surface MSD 
we find at short times (t t K ) 



(z 2 (t)} - 2D s t 



2 D h 



3y^D s \t 



t 



1/2 



{*\t)) n . 



(12) 

In this regime, that is, the diffusion to leading order is 
confined to the surface and exchange with the bulk leads 
to a higher order correction, when D s <C Db, that invokes 
~ t c 'I 2 superdiffusion (see below). Conversely, for longer 
times f K <Ct<f„, N s (t) decays perceptibly and 



{t\t))n 



2Dd h 



--D h Jtt K 



2yfKD s y r tt^+2Dbt 



(13a) 
(13b) 



At even longer times t a -C t -C fa we see the influence of 
the cylindric geometry in the logarithmic dependencies 



(z\t)) 



°^]og(JL) (14a.) 



2D b t (14b) 




where logC = 7 = 0.5772 is Euler's constant. In this 
regime (z 2 {t)) thus exhibits a saturation unique to the 
cylinder case. Finally, at very long times t 3> fa due to 
radial equilibration a linear diffusive behavior yields 



(z 2 (t)) r 



^D s t + 4D b vM* t 



4£>, 



tb 



fa 

-t + 2D b t, 



(15a) 
(15b) 
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Figure 2: Mean squared displacement showing the various 
effective surface diffusion regimes. Note the transient plateau. 



i.e., the combination of bulk and surface diffusion gives 
rise to an effective diffusivity involving all time scales. 

The behavior of the surface MSD (z 2 (t)) in absence of 
surface diffusion (D s = 0) is shown in Fig. [2] obtained 
from numerical Laplace inversion of Eqs. (fTTj) . ©, and 
©. Note that the time scales (in dimensionless units) 
t R = 1CT 6 , t a = 1, and tb = 10 6 were chosen to be 
well separated, to visualize the four scaling regimes. The 
other parameters were a = 5, Db = a 2 /t a = 25, k = 
aj\Jt a t K = 5 x 10 3 , and b = a^/t b /t a = 5 x 10 3 . 

Surface propagator. Our formalism also produces the 
exact distributions C(r,t) and n(z,t). For times longer 
than tb we already found that the diffusion behavior is 
normal albeit with rescaled diffusivity and the surface 
propagator turns into a Gaussian. We recover a Gaussian 
for n(z, t) in absence of bulk coupling, n = 0. For times 
shorter than tb and k > we observe deviations from the 
Gaussian form. In most relevant cases D s <C Db and we 
therefore consider the case D s = in what follows. 

At shorter times t <C t a < t b we have ut a 3> 1 and 
utb ^> 1 and qa — a^Jk 2 + u/Db > a^J u/Db = y/ut a 

A and 



a\Jk 2 + u/Db > a.\J u/Db = 
1 such that qa ^> 1 and qb 3> 1. Therefore Ai ? 
the propagator reduces to 



n(k, u) 



N 



+ Ky/k 2 + U/D b 



(16) 



In the central part of the DF where k 2 > u/Db (i.e., z < 
yfDbt), we obtain the normalized Cauchy distribution 



n(k, u) 



u + n\k\ 



<^ n(z,t) 



N Kt 



ir (z 2 + kH 2 ) ' 



(17) 



This interesting result is analogous to the findings from 
Ref. [2| for a flat surface found from involved scaling argu- 
ments. In contrast, here we derive the Cauchy law from 
an exact approach allowing us to study the transition to 
other regimes explicitly. To do so we introduce the range 



of validity £c(t) — \/Dbt of the Cauchy region. While 
at distances z > £c(t) we observe a Gaussian cutoff, for 
z < lc{t) the Cauchy approximation is valid. From this 
Cauchy part we obtain the superdiffusive contribution 



2 nt dz 



la® ^ 

-l C (t)K(z 2 + K 2 t 2 ) 



2nJD~ h t i/2 



(18) 



to the MSD, consistent with the exact result (TT2"]) [with 
D s = 0]. Calculation of the MSD from the full Eq. © 
however requires the k — » limit and thus the extreme 
wings of the distribution are explored. As the system 
evolves in time the central Cauchy part spreads. As al- 
ready in the regime t K < t < t a we have Dbt < K 2 t 2 the 



asymptotic behavior 



can no longer be observed. 



At intermediate times t K <C t <C t a we find n(k, u) 
Nq sJt K /{u + Dbk 2 ) corresponding to the Gaussian 



n(z,t) = N a (t) 



1 



4nD b t 



exp 



AD h t) 



(19) 



with N s (t) = f_ oo n(z,t)dz = N oy /t K /{-Kt). For longer 
times t a <C t tb the propagator acquires the shape 
n(k, u) = -aN log (C 2 [a 2 k 2 + ut a }/4) /[2k] which yields 
Eq. (JTSJ with N s (t) = N iftJ^/[2t}. Finally, for very 
long times t t b we again find Eq. lfl9|) with the satu- 
ration value N s (t) = 2No^/t a t K /tb. Indeed one can show 
by the shift theorem of the Laplace transformation that 
any function of the argument u + Dbk 2 will lead to the 
Gaussian shape (|19[) . with appropriate normalization. 

First passage. The DF of times p(t) a particle spends 
in the bulk after detachment from the cylinder can be 
calculated explicitly (details of the calculation will be 
presented elsewhere [12|). To this end we initially place 
the test particle at radius a < ro < b and calculate when 
it is first adsorbed at r = a. With to = r\/Db this first 
passage problem defines p(t) by 



p(t) = 2na 



D, 



dC(r,z,t) 



Or 



dz 



(20) 



in terms of the radial flux into r 
form of p(t) then becomes 



p{u) 



a. The Laplace trans- 

h {Vutb) Kg (y/uto~) + Ki (yfutl) Ig (y/uto~) 

h (Vutb) K Q (s/uQj + Ki (s/vfo) J (\/uQ) 



— — (21) 

For tb — > oo, p(u) ~ Ko (\/uto) / K (\/ut a ) while for 
ro = a we recover the sharp form p(t) = S(t) as it should. 
At shorter times t -C t a < to we obtain the expansion 



P(t) 



'Y) 



(r<l- 



1 + 



D h t 



(22) 



r V^D b t 3 V 4ar 

This is to leading order the first passage DF for a one- 
dimensional random walk reweighted by the ratio y/ a / r o- 
Keeping the distance ro — a fixed but letting both ro and 
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a tend to infinity we recover the result for a flat surface 
for which the ID first passage remains valid at all times. 

At longer times t a <C t <C % the logarithmic form 
p(u) ~ 1 — 2 logfrn la) / log(l / \ut a ] ) yields. From Taube- 
rian theorems [l3j we infer the first passage behavior into 
a cylinder of radius a, 



Pit) 



21og(r /a) 
t\og 2 (t/t a y 



(23) 



see also Refs. [ll|, 
return times to the 
fusive coupling z 2 : 
to the DF X(z) ~ l/(zlog 2 z) of effective dislocations 
along the cylinder, 
the outer cylinder 
exponential cutoff, 



14j j . We note that a distribution of 
cylinder of this form implies by a dif- 
i t that a single bulk excursion leads 



Finally in the long time limit t i& 
comes into play and p(t) attains an 
leading to a mean first passage time 



(t) =b 2 log{r /a)/[2D b }. 



(24) 



Extremely long jumps. An interesting behavior occurs 
when we relax the strong coupling condition t K <C t a . 
We said that when t ^S> t a the diffusing particles feels a 
change in geometry from planar to cylindrical, and the 
time scale t c = y/t a t K = a/n occurs in our expressions. 
If we have t a < t < t c < t b then N s (t) « A and 



(z 2 (i)) =2D s t + 



2D h t 2 



t c \og 2 [U/iCHa)] 1 



(25) 



a ballistic behavior with logarithmic correction: For 
D s <C -D& the superdiffusion is even stronger than for 
t <C t a . Similarly, the propagator reduces to 



n(k, u) = N Q /[u + nqKi(aq) / K (aq)]. 



(26) 



In the range a <C z <C \JD b t (i.e., ^/ut a <§; afc <C 1), 
KqKi(aq)/K (aq) ~ (/t/a) log[2/(Cafe)], an extremely 
weak fc dependence: this pre-cutoff tail of n(fc, u) is heav- 
ier than any normalizable power law, due to the heavy 
tailed distribution of bulk mediated dislocations. 

Discussion. We established an exact approach to 
BMSD on a reactive cylindrical surface revealing four 
distinct surface diffusion regimes. In particular our for- 
malism provides a stringent derivation of the transient 
superdiffusion discussed earlier and explicitly quantifies 
the transition to other regimes. Notably we revealed a 
saturation regime for the surface MSD that becomes rel- 
evant at times above which the diffusing particle feels 
the curvature of the cylinder surface (t a ). This behavior, 
caused by the cylindrical geometry, stems from an in- 
teresting balance between a net flux of particles into the 
bulk and the fact that particles with a longer return time 
also lead to an increased effective surface relocation. In 
absence of an outer cylinder the saturation is terminal, 
while in its presence the surface MSD returns to a linear 
growth in time. This observation will be important in fu- 
ture models of BMSD around cylinders and particularly 



for the interpretation of experimental data obtained for 
BMSD systems. We note that in the proper limit a — > oo 
the previous results for a planar surface are recovered. 
Relaxing the strong coupling condition we demonstrated 
the existence of an almost ballistic surface diffusion be- 
havior, a case that might be relevant for transport along 
thin cylinders such as DNA. 



In Ref. [ll| it was shown that the scaling behavior in 
the regimes below and above t a can be probed experi- 
mentally by NMR methods measuring the BMSD of wa- 
ter molecules along imogolite nanorods over three orders 
of magnitude in frequency space. For larger molecules 
such as a protein of approximate diameter 5 nm we ob- 
serve a diffusivity of 10 _6 cm 2 /sec such that for instance 
the saturation plateau around a bacillus cell (radius 1/2 
/xm) sets in at around t a = 2.5 msec which might give 
rise to interesting consequences for the material exchange 
around such cells. In general, the relevance of the indi- 
vidual regimes will crucially depend on the scales of the 
surface radius and the diffusing particle (and therefore 
its diffusivity). It was discussed previously that even the 
superdiffusive short-term behavior may become relevant 
[H d, EH- In general, in a given system the separation 
between the various scaling regimes may not be sharp. 
Moreover typically a single experimental technique will 
not be able to probe all regimes. It is therefore vital to 
have available a solution for the entire BMSD problem. 
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